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Abstract 

We analyze the hitting time distributions of stock price returns in different time windows, char- 
acterized by different levels of noise present in the market. The study has been performed on two 
sets of data from US markets. The first one is composed by daily price of 1071 stocks trade for the 
12-year period 1987-1998, the second one is composed by high frequency data for 100 stocks for 
the 4-year period 1995-1998. We compare the probability distribution obtained by our empirical 
analysis with those obtained from different models for stock market evolution. Specifically by fo- 
cusing on the statistical properties of the hitting times to reach a barrier or a given threshold, we 
compare the probability density function (PDF) of three models, namely the geometric Brownian 
motion, the GARCH model and the Heston model with that obtained from real market data. We 
will present also some results of a generalized Heston model. 
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I. INTRODUCTION 

The interest of physicists in intedisciplinary researches has been largely increased in re- 
cent years and one of the developing field in this context is econophysics. It applies and 
proposes ideas, methods and models in statistical physics and physics of complex systems to 
analyze data coming from economical phenomena jjj. Several statistical properties verified 
in financial quantities such as relative price changes or returns and their standard deviation, 
have enabled the establishment of new models which characterize systems ever better j^J. 
Moreover the formalism used by physicists to analyze and to model complex systems consti- 
tutes a specific contribution that physics gives to many other fields. Complex systems in fact 
provide a very good paradigm for all those systems, physical and non-physical ones, whose 
dynamics is driven by the nonlinear interaction of many agents in the presence of " natural" 
randomness [3] . The simplest universal feature of financial time series, discovered by Bache- 
lier j^J, is the linear growth of the variance of the return fluctuations with time scale, by 
considering the relative price changes uncorrelated. The availability of high frequency data 
and deeper statistical analyses invalidated this first approximated model Q], which is not 
adequate to catch also other important statistical peculiarities of financial markets, namely: 
(i) the non-Gaussian distribution of returns, (ii) the intermittent and correlated nature of 
return amplitudes, and (iii) the multifractal scaling 5j, that is the anomalous scaling of 
higher moments of price changes with time. 

In this paper we focus our attention on the statistical properties of the first hitting 
time (FHT), which refers to the time to achieve a given fixed return or a given threshold 
for prices. Theoretical and empirical investigations have been done recently on the mean 
exit time (MET) |f| and on the waiting times [jj] of financial time series. We use also the 
term " escape time" to include the analysis of times between different dynamical regimes 
in financial markets done in a generalized Heston model jsj. Markets indeed present days 
of normal activity and extreme days where high price variations are observed, like crash 
days. To describe these events, a nonlinear Langevin market model has been proposed in 
Refs. H, where different regimes are modelled by means of an effective metastable potential 
for price returns with a potential barrier. We will discuss three different market models 
evidencing their limits and features, by comparing the PDF of hitting times of these models 
with those obtained from real financial time series. Moreover we will present some recent 
results obtained using a generalized Heston model. 
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II. MODELS FOR STOCK MARKET EVOLUTION 
A. The geometric random walk 

The most widespread and simple market model is that proposed by Black and Scholes to 
address quantitatively the problem of option pricing Q]. The model assumes that the price 
obeys the following multiplicative stochastic differential equation 

dp(t) = fi ■ p(t) ■ dt + a • p(t) ■ dW{t) (1) 

where /i and a are the expected average growth for the price and the expected noise intensity 
(the volatility) in the market dynamics respectively, dp jp is usually called price return. The 
model is a geometric random walk with drift fi and diffusion a. By applying the Ito's lemma 
we obtain for the logarithm of the price 
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d lnp{t) = (fj, - y) • dt + a ■ dW{t) . (2) 

This model catches one of the more important stylized facts of financial markets, that is the 
short range correlation for price returns. This characteristic is necessary in order to warrant 
market efficiency Short range correlation indeed yields unpredictability in the price time 
series and makes it difficult to set up arbitrage strategies. 

The statistical properties of escape times r for this process are well known, and the PDF 
of escape time r, F(r,p ), was obtained analytically 11]. If the starting value of the price 
is po at time t = 0, the distribution of the time r to reach a barrier at position h is given by 
the so called inverse Gaussian 



F{r,p ) = t P ° i • ex P 
V ZTia z T A 



(h-po- fir) 
2a 2 r 



(3) 



which is well known among finance practitioners to price exotic options like barrier op- 
tions or to evaluate the probability for a firm to reach the zero value where it will 
remain forever. The shape of the distribution is shown in Fig. for two different cases. The 
asymptotic expressions of PDF show in one case a power law tail with exponent -1.5 
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FIG. 1: Inverse Gaussian distribution obtained with the parameters a = 1.0, {h — po) = 1.0, for 
/i = (solid line) and /i = 1 (dotted line). Inset: log-log plot of the same PDF. 

The distribution of hitting times for price returns instead is very simple. In the geometric 
Brownian motion the returns are independent, so the probability to observe a value after a 
certain barrier is given by the probability that the "particle" doesn't escape after n — 1 time 
steps, multiplied by the escape probability at the n th step 

F r (r) = (l-p)-p n ' 1 = ((l-p)-exp[(n-l)\np}, n=^ (6) 

where p is the probability to observe a return inside the region defined by the barrier, At 
is the observation time step and r is the escape time. So the probability is exponential 
in time. The geometric Brownian motion is not adequate to describe financial markets 
behavior, because the volatility is considered as a constant parameter and the PDF of the 
price is a log-normal distribution. 



B. The GARCH and the Heston models 



Price returns have indeed properties that cannot be reproduced by the previous simple 
model: (i) price return distribution has fat tails; (ii) price returns have short range correlation 
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but the volatility is a stochastic process with long range correlation 
variability in time of the volatility indeed depends not only on the Fundamentals of the 
firm but also on the market conditions. Volatility is usually higher during crisis periods 
and has also an almost deterministic intra-day pattern during the trading day, being higher 



near market opening and closure. So the volatility can be considered as a stochastic process 
itself and it is characterized by long range memory and clustering. More realistic models to 
reproduce the dynamics of the volatility have been developed. Here we will present two of 
them: the GARCH and the Heston models. 

The GARCH(p,q) process (generalized autoregressive conditional heteroskedasticity), 
which is essentially a random multiplicative process, is the generalization of the ARCH 
process and combines linearly the actual return with p previous values of the variance and 
q previous values of the square return [RJ. The process is described by the equations 

of = a + a x x\_ x H V a 2 q x 2 _ q + fao 2 t _ x H H Pp°t- P > x t = Vt ' °u ( 7 ) 

where a.i and fa are parameters that can be estimated by means of a best fit of real market 
data, Xt is a stochastic process representing price returns and is characterized by a standard 
deviation o t . The GARCH process has a non-constant conditional variance but the variance 
observed on long time period, called unconditional variance, is instead constant and can 
be calculated as a function of the model parameters. It has been demonstrated that x\ of 

re autoco- 
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GARCH(1,1) is a Markovian process with exponential autocorrelation, while t 
variance of GARCH(p,q) model is a linear combination of exponential functions 
will consider the simpler GARCH(1,1) model 



of = a + {atirfi^ + fa)<J 2 -ii x t = r] t ■ cr t . (8) 

The autocorrelation function of the process x t is proportional to a delta function, while 
the process xf has a correlation characteristic time equal to r =| ln{a,\ + fa) | _1 and the 
unconditional variance equal to a 2 = cno/(l — a\ — fa). So it is possible to fit the empirical 
values of these two quantities by adjusting few parameters. Specifically a\ and fa regulate 
the characteristic time scale in the correlation, while can be adjusted independently to 
fit the observed unconditional variance. 

In the Heston model Q| the dynamics is described by a geometric Brownian motion 
coupled to a second stochastic process for the variable v = a 2 . The model equations are 

dlnp{t) = (p -'!&-). dt + y/v(fij ■ dW(t), (9) 
dv{t) = a(b-v(t))-dt + c^/v(t)-dZ(t), 

where W(t) and Z(t) are uncorrelated Wiener processes with the usual statistical properties 
< dWi > = 0, < dWi(t)dWj(t') > = 8{t—t') 5ij, but can be also correlated The process 



for v is called Cox-Ingersoll-Ross JLfJ and has the following features: (i) the deterministic 
solution tends exponentially to the level b at a rate a (mean reverting process); (ii) the 
autocorrelation is exponential with time scale r = a -1 . Here c is the amplitude of volatility 
fluctuations often called the volatility of volatility. Once again we have a model with short 
range correlation that mimics the effective long range correlation of the markets using large 
values of r. The process is multiplicative and values of v can be amplified in few steps 
producing bursts of volatility. If the characteristic time is large enough, many steps will be 
required to revert the process to the mean level b. So the longer the memory is, the longer 
the burst will survive. For little correlation times the process fluctuates uniformly around 
the mean level b, whereas for large correlation times v presents an intermittent behavior with 
alternating activity of burst and calm periods. The model has been recently investigated by 
econophysicists lj| and solved analytically 

The two models presented so far are a more realistic representation of financial market 
than the simple geometric Brownian motion, even if they do not reproduce quantitatively 
the form of the long time correlation observed for the volatility. We use a set of 1071 daily 
stock price returns for the 12-year period 1987-1998, and we compare the results obtained by 
simulation of the GARCH and Heston models with those obtained from real market data. 
The parameters in the models were chosen by means of a best fit, in order to reproduce 
the correlation properties and the variance appropriate for real market. Specifically for the 
GARCH model we used values a± = 0.07906 and Pi = 0.90501 obtained elsewhere [l^ to 
fit the correlation time of daily market data, and «o = 7.7 • 10 -6 in order to fit the average 
standard deviation of our data using the formula for unconditional variance presented in 
the previous section. For the Heston model we used a = 4.5 • 10~ 2 , b = 8.62 • 10~ 5 and 
H = 5.67 ■ 10~ 4 , obtained in a recent work suitable for daily returns and c = 10.3 ■ 10~ 3 , 
as before, to fit the average standard deviation of our data set. Using these parameters we 
obtain distribution for price returns that are in reasonable agreement with real market data, 
as shown in Fig. 121 The two models approximate the return distributions of real data better 
than the Gaussian curve. In particular the Heston model gives the best agreement. The 
chosen parameter set therefore is good enough to fit the dynamics of our data. 

In order to investigate the statistical properties of escape times r we choose two thresholds 
to define the start and the end for the random walk. Specifically we calculate the standard 
deviation <jj, with i — 1, . . . , 1071 for each stock over the whole 12- year period. Then we set 
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FIG. 2: Probability density function of stock price returns for: real market data (circle), GARCH 
model (square) and Heston model (diamond). The black solid line is a Gaussian distribution having 
the same standard deviation of real data. 

the initial threshold to the value 0.1-Cj and as final threshold the value — 2-crj. The thresholds 
are different for each stock, the final threshold is considered as an absorbing barrier. The 
resulting experimental distribution reported in Fig. El has an exponential tail but it deviates 
from the exponential behavior in the region of low escape times. Specifically low escape times 
have probability higher than the exponential. We recall that for the geometric Brownian 
motion model this distribution should be exponential over the entire r axis. So the first 
conclusion we can draw from our analysis is that the basic geometric Brownian motion is 
not adequate to explain the distribution of r. 

In order to reproduce more closely the situation present in real market we choose a s t ar t 
only once, specifically we place the random walker in the initial starting position and we 
set the initial volatility value. When the random walker hits the barrier, we register the 
time and we place the walker again in the initial position, using the volatility of the barrier 
hitting time. So the random walker can experience different initial volatility values as in 
real markets. The results are reported in the two panels of Fig. Olfor both the GARCH and 
the Heston models, and we see that these models provide a better agreement with real data 
than the geometric Brownian motion. Moreover for the GARCH model the agreement is only 
qualitative, whereas the Heston model is able to fit the empirical distribution quantitatively. 
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FIG. 3: PDF of escape times of price returns for the (a) GARCH and (b) Heston models (diamond) 
compared with the distribution obtained from real market data (circle). The process starts at 
(— 0.1(7 ) = —0.0022, the barrier is located at (—2. Oct). The straight line is an exponential fit of the 
distribution tail. 

C. The modified Heston model 



Here we consider a generalization of the Heston model, by considering a cubic nonlinear- 
ity. This generalization represents a fictitious " Brownian particle" moving in an effective 
potential with a metastable state, in order to model those systems with two different dynam- 
ical regimes like financial markets in normal activity and extreme days |9]. The equations 
of the new model are 

dx(t) = + ^ dt + yfitf) dW 1 (t) (10) 

dv(t) = a(b-v(t)) dt + cy/v(£) dW 2 {t), (11) 

where U(x) = 2x 3 + 3x 2 is the effective cubic potential with a metastable state at x me = 0, a 
maximum at xm = — 1, and a cross point between the potential and the x axes at xj = — 1.5. 
In systems with a metastable state like this, the noise can originate the noise enhanced 
stability (NES) phenomenon, an interesting effect that increa ses, instead of decreasing, the 



stability by enhancing the lifetime of the metastable state [19|,|20j. The mean escape time r 
for a Brownian particle moving throughout a barrier AU, with a noise intensity v, is given 
by the the well known exponential Kramers law r = exp [AU /v], where r is a monotonically 
decreasing function of the noise intensity v. This is true only if the random walk starts 
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from initial positions inside the potential well. When the starting position is chosen in 
the instability region x D < xm, t exhibits an enhancement behavior, with respect to the 
deterministic escape time, as a function of v. This is the NES effect and it can be explained 
considering the barrier "seen" by the Brownian particle starting at the initial position Xq, 
that is AUin = U(x max ) — U(xq). In fact AUi n is smaller than AU as long as the starting 
position xq lies in the interval / = [x/,xa/]. Therefore for a Brownian particle starting 
from an unstable initial position, from a probabilistic point of view, it is easier to enter into 
the well than to escape from, once the particle is entered. So a small amount of noise can 
increase the lifetime of the metastable state. For a detailed discussion on this point and 
different dynamical regimes see Refs. J^oJ- When the noise intensity v is much greater than 
AU, the Kramers behavior is recovered. 

Here, by considering the modified Heston model, characterized by a stochastic volatility 
and a nonlinear Langevin equation for the returns, we study the mean escape time as a 
function of the model parameters a, b and c. In particular we investigate whether it is 
possible to observe some kind of nonmonotonic behavior such that observed for r vs. v in 
the NES effect with constant volatility v. We call the enhancement of the mean escape 
time (MET) r, with a nonmonotonic behavior as a function of the model parameters, NES 
effect in the broad sense. Our modified Heston model has two limit regimes, corresponding 
to the cases a = 0, with only the noise term in the equation for the volatility v(t), and 
c = with only the reverting term in the same equation. This last case corresponds to the 
usual parametric constant volatility regime. In fact, apart from an exponential transient, 
the volatility reaches the asymptotic value b, and the NES effect is observable as a function 
of b. To this purpose we perform simulations by integrating numerically the equations (10) 
and (jllj) using a time step At = 0.01. The simulations were performed placing the walker 
in the initial positions x located in the unstable region [xj,xm] and using an absorbing 
barrier at x = —6.0. When the walker hits the barrier, the escape time is registered and 
another simulation is started, placing the walker at the same starting position xq, but using 
the volatility value of the barrier hitting time. 

The mean escape time as a function of b is plotted in Fig. 0] for the 3 different starting 
unstable positions and for c = 0. The curves are averaged over 10 5 escape events. The non- 
monotonic behavior is present. After the maximum, when the values of b are much greater 
than the potential barrier height, the Kramers behavior is recovered. The nonmonotonic 
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FIG. 4: Mean escape time r for 3 different unstable starting positions, when only the reverting 
term is present: a = 10~ 2 , c = 0). 



behavior is more evident for starting positions near the maximum of the potential. For a = 
the system is too noisy and the NES effect is not observable as a function of parameter c. 
The presence of the reverting term therefore affects the behavior of r in the domain of the 
noise term of the volatility and it regulates the transition from nonmonotonic to monotonic 
regimes of MET. The results of our simulations show that the NES effect can be observed 
as a function of the volatility reverting level b, the effect being modulated by the parameter 
(ab)/c. The phenomenon disappears if the noise term is predominant in comparison with 
the reverting term. Moreover the effect is no more observable if the parameter c pushes the 
system towards a too noisy region. When the noise term is coupled to the reverting term, we 
observe the NES effect on the variable c. The effect disappears if b is so high as to saturate 
the system. 

We compare now the theoretical PDF for the escape time of the returns with that obtained 
from the same real market data used in the previous section. We define two thresholds, 
Axi = O.Iuax and Axj = 1.0<j& x , which represent respectively start point and end point 
for calculating r. The standard deviation a^ x of the return series is calculated over a 
long time period corresponding to that of real data. The initial position is xq = —1.25 
and the absorbing barrier is at x a b s = —6.0. For the CIR stochastic process v we choose 
Vstart = 8.62 • 10 -5 , a = 10 _1 , b = 4.5 and c = 2 • 10 _1 . The agreement with real data is very 
good. At high escape times the statistical accuracy is worse because of few data with high 
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FIG. 5: Probability density function of the escape time of the returns from simulation (solid line), 
and from real data (black circle). 

values. The parameter values of the CIR process for which we obtain this good agreement 
are in the range in which we observe the nonmonotonic behavior of MET. This means that 
in this parameter region we observe a stabilizing effect of the noise in time windows of prices 
data for which we have a fixed variation of returns between Ax, and Axf. This encourages 
us to extend our analysis to large amounts of financial data and to explore other parameter 
regions of the model. 




III. ESCAPE TIMES FOR INTRA-DAY RETURNS 



In this last section we discuss the results obtained with the same analysis of the previous 
section, applied to a different data set at intra-day time scale. The data set contains 100 
stocks in the 4-year period 1995-1998. The stocks considered are those used, in that period, 
in the S&P100 basket. We are dealing therefore with highly capitalized firms. The data 
are extracted from the Trade and Quote database. The stocks are distributed in different 



market sectors as illustrated in Ref. 



21(. For the analysis we considered the return on a 



time interval equal to St = 1170sec, which is approximately equal to 20 minute and it is 
contained in a market day exactly 20 times. So we have 19 price returns per day, which 
amounts to 20220 points in the whole period of 4 years per each of the 100 stocks. We 
used the value —0.5 • cr as a start position and the value —7.0 • o as absorbing barrier. We 
can choose a so high value for the barrier because return distribution on intra-day time 
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scale have tails fatter than daily return distribution, therefore the statistical accuracy for so 
high barrier value is good enough. The distribution of escape times obtained is reported in 
Fig. |Ulin a semi-logarithmic plot. It has an exponential trend superimposed to a fluctuating 
component. One can recognize that the period of the fluctuation is 1 trading day, so this 
effect has to be ascribed to something that happens inside the daily activity. To describe 
better this aspect we record, for each barrier hitting event, the hour when the event occurs 
and we build a histogram showing the number of barrier hitting events as a function of the 
day time. This histogram (see the inset of Fig. IHJ) clearly shows that the barrier hitting takes 

Itf 3 

gio" 4 

time (min) 

FIG. 6: PDF of escape times of price returns for intra-day price returns obtained from real market 
data for the 100 stocks of the S&P100 basket in the 4-year period 1995-1998. The process starts at 
—0.5<t, the barrier is located at —7.0a. Inset: Frequency distribution of the barrier hitting event, 
within a day, as a function of the hour. The internal inset shows the volatility observed in the 
same moment of the trading day. 

place more frequently near the opening and the closure of the market. This happens because 
the volatility follows a well known deterministic pattern during the day, being higher near 
the market opening and closure, and lower in the middle of the trading day. In the same 
inset we report an estimation of the volatility per hour, which we calculate as the standard 
deviation of the return observed in that hour, in the whole period for all the 100 stocks. 
The figure shows that the volatility has a pattern reproducing that observed in the barrier 
hitting event histogram. This is in agreement with our considerations. 
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IV. CONCLUSIONS 



We studied the statistical properties of the hitting times in different models for stock 
market evolution. We discussed limitations and features of the basic geometric Brownian 
motion in comparison with more realistic market models, such as those developed with a 
stochastic volatility. Our results indeed show that to fit well the escape time distribution 
obtained from market data, it is necessary to take into account the behavior of market 
volatility. In the generalized Heston model the reverting rate a can be used to modulate 
the intensity of the stabilizing effect of the noise observed (NES), by varying b and c. In 
this parameter region the probability density function of the escape times of the returns fits 
very well that obtained from the real market data. The analysis on intra-day time scale 
shows another peculiarity: the intra-day volatility pattern produces periodic oscillations in 
the escape time distribution. This characteristic will be subject of further investigation. 
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